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Abstract. - We consider a tunnel junction between two arbitrary non-linear systems in any 
dimension, which can be different. We show that the tunneling charge can be detected using 
three alternative methods based on current measurements. Besides being technically easier com- 
pared to noise measurements, these methods present valuable advantages: they do not require 
the knowledge of the underlying models, and some are accessible in the experimentally convenient 
low-voltage regime, where heating effects are reduced. The first method is based on the AC con- 
ductance, while the two others are based on photo-assisted current (PAC) and can be implemented 
for any time-dependence of the tunneling amplitude. These are promising for edge states in the 
regime of the fractional quantum Hall effect (FQHE): the Hamiltonian does not have to be speci- 
fied and can incorporate non-universal interactions between the edges, and it is more convenient 
to use an AC gate voltage rather than an AC bias. These methods apply for instance to weak 
barriers in 1-D systems, Superconductor-Insulator-Normal (SIN) or graphene-like structures. 



One of the most exciting hallmarks of electronic interac- 
tions are elementary excitations with non-integer charges 
or with a double charge in the super-conducting state 
[l][2]. Searching for methods to measure them has at- 
tracted tremendous efforts on both the theoretical and 
experimental level. Fractional charges have even been pre- 
dicted recently in graphene-like structures [s] or topolog- 
ical insulators [4|. A rich variety has been mainly offered 
by the FQHE, where quasi-particles are predicted also to 
obey fractional statistics. Their measurement plays a cru- 
cial role as it may shed light on the underlying physics. 
While for simple filling factors (defined as the dimension- 
less ratio between the number of electrons and that of the 
flux quanta) v = 1/(1 + 2n), with n integer, the charge is 
simply q = ve, it was claimed that for the Jain series [s], 
u = p/{l -\- '2np), the fundamental charge is e/(l -I- 2np). 
Recently, there has been a resurgence of interest in the 
non-abelian = 5/2 fractional state, whose different de- 
scriptions have a common prediction: the existence of a 
charge e/4. Its detection, being one of the crucial steps 
toward revealing non-abelian statistics [6], has not been 
conclusive yet [7||9]. Indeed, for any filling factor, theory 
and experiment have not been fully reconciled so far. 

Theoretical suggestions to determine the charge of car- 
riers has been mostly based on the Poissonian formula for 



the noise S at sufficiently weak current /, 



S = ql, 



(1) 



obtained for some specific systems 10 



the charge of the carriers experimentally II 



and used to access 
It turns out 



that this formula does not depend on the choice of the 
underlying model 12 13 . However, its experimental ap- 
plications encounter a number of difficulties. First, noise 
is more difficult to measure compared to current. Second, 
Eq. ([!]) requires the condition of high voltage compared 
to temperature, qV ^ T (we set fc^ = 1 and h — 1) 
[14] , which can cause heating effects in experiments. In 
Heiblum's group 



15 



where the regime of higher volt- 
ages is investigated, this equation has been replaced by 
a conjectured one: that for non-interacting particles with 
charge q. This has led to results which have not been fully 
understood, even for simple fractions. There has also been 
a theoretical proposal to use the photo-assisted noise to de- 
tect the charge at simple fractions 16 , a result which we 



can indeed generalize to any filling factors and arbitrary 
systems in any dimension 13 . But still such a proposal 



is based on noise and requires qV » T. 

On the other side, there are some experiments based 
on conductance measurements [sj. Their disadvantage is 
that they rely on a particular model for the edges, thus on 
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factors such as the edge reconstruction and Coulomb inter- 
actions. For instance, the exact structure of edge modes 
in "hole conjugate" sates, at fillings 2/3, 3/5, etc., remains 
still an unsolved issue pT?, . 

Here we propose alternatives to measure the charge 
which avoid one or, simultaneously, many of these limita- 
tions: they are based on current measurements, and some 
of them are possible in the low voltage regime, which re- 
duces heating. The important aspect of these results is 
their universality, i.e., the independence of the model as 
long as the system is kept in a weak tunneling regime. We 
exclude only the treatment of a super-current, which in an 
Superconductor-Insulator-Superconductor (SIS) structure 
can be suppressed by a magnetic field. More generally, 
we can allow for two different systems, for an STM and 
another system to probe, and for coupling to any type of 
electromagnetic environment. Our computation applies to 
both tunneling and weak backscattering in correlated 1- 
D systems. The charge q enters the equations (|8 17) in 
many different ways, which allows us to propose a variety 
of methods to detect it. Some of them require zero-bias 
anomaly (ZBA), such as in 1-D systems and the FQHE, 
systems with disorder and interactions and/or coupled to 
an ohmic environment. 

Two of our proposals are based on the Photo- Assisted 
Current (PAC) 



16 18-20 



Indeed, we make another 
important achievement: it gives a generalization of the 
theory of Tien-Gordon jl^ intended to deal with Photo- 
Assisted Tunneling (PAT) in SIS junctions, and that by 
Tucker 21 . We extend these theories to treat: 



• A general Hamiltonian required only to conserve the 
charges in both systems (in the absence of tunneling), 
thus allowing for Coulomb interactions between them. 
This is of great importance in the FQHE. 

• A tunneling charge different from e, without specifying 
the explicit form of the tunneling term. 

• Differentials of the current with respect to the Fourier 
components of an "effective" voltage kept finite. 

• Arbitrary periodic dependence on time of the voltage 
or/and the tunneling amplitudes. 

The latter is easier to settle in the FQHE. Indeed, time- 
dependent barriers in one dimension have been studied 




Fig. 1: Tunnel junction between two systems 1 and 2 in any 
dimension with arbitrary interactions or disorder, subject to 
time-dependent potential Vi,2(i). They can be either similar or 
different, such as fractional edge states with equal or different 
filling factors, or SIN structures. One can also view one system 
as a probe, such as STM tip or a Fermi Hquid coupled to edge 
states or to a quantum dot. Mutual Coulomb interactions are 
allowed, we require only Qi and Q2 to commute with the total 
Hamiltonian 'Ho in the absence of tunneling. 



tunneling between the two systems by 'Ht- Thus the total 
Hamiltonian reads : 



Q 



(2) 



We allow for an arbitrary dependence of the tunneling 
amplitudes Tk-i^k2i't) on the states ki and fc2 on both sides, 
as well as for their time-dependence, assuming it only to 
be identical for the different states, i.e., the dimensionless 
function: 7(f) = r/ci,fe2(^)/rfci,fc2(0) does not depend on 
fci,fc2. This is reasonable if time dependence were to be 
controlled by a nearby gate. Therefore, denoting by T the 
tunneling operator at time 0, we have: 



Hr-7(^)r + 7*(^)r^ 



(3) 



where we do not have to write the explicit form of T. For 
the perturbative computation carried out later, we require 



previously 23 , where the current was expressed explicitly, only that (T(t)T(O)) = in the absence of tunneling 24 



depending on the underlying model, thus have a different 
spirit and aim. Besides, we offer a unified background 
which contains all those different results, and make the 
bridge with the PAT. 

Model: We start from two (identical or different) sys- 
tems whose Hamiltonian Hq incorporates interactions or 
disorder which do not need to be specified. In particular, 
we do not require Ho to be decoupled into two Hamiltoni- 
ans, thus can include mutual Coulomb interactions. The 
only requirement for Jio is that it does not change the total 
charges Qi and Q2- (see Fig. [T]) [Q/,Ho] = for / = 1,2. 
We describe a gauge invariant coupling to the applied po- 
tentials Vi{t) and V2it) by Hq = QiVi{t) + (52^^2(0, and 



which implies that the super-current is zero, or has to be 
negligible compared to that of quasi-particles. In an SIN 
structure, one has q = 2e (respectively, q — e) at low 
(respectively, high) energies compared to the gap. In the 
FQHE as well, q can depend on energy scales. 

The operator T transfers a charge q, thus translates the 
charge difference Q = {Q2 — Qi)/2 hy q. Therefore it con- 
tains, without any loss of generality, an exponential of a 
displacement operator conjugate to Q: T — T'e''*, where 
[$, Q] = i, and T' commutes with both $ and Q. This is in 
particular suited to include coupling to an electromagnetic 
environment or even to another quantum system. Using 
the scalar phase (p defined as </> = Vi{t) — V2{t) = V{t), 
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and the standard time-dependent unitary transformation 
U{t) = e-'O-^^*), one has U{t)TU'^{t) = e'^-^^'^T. This al- 
lows us to absorb Hq by applying V. — UHW — iUdthl^ = 
Hq + Hrit), so that only Ht, being gauge non-invariant, 
acquires additional time-dependence: 



7(t)|. 



(4) 
(5) 



The phase (j>{t) and V emanate from both V{t) and 
arg7(i). More precisely, the acquired phase arg['y{t)]/q + 
4){t) = Vt-\- (t)(t) is decomposed so that its time derivative, 
yielding an "effective" voltage, contains a DC part V and 
a time-dependent one V{t) = dt(l>{t). 

Now it is only 'Hj'{t) which drives the system from an 
equilibrium state. Treating it as a perturbation, we apply 
the standard procedure of expanding the evolution oper- 
ator in the interaction representation, keeping only the 
first non-vanishing contribution to /, the average of the 
tunneling current operator: 



/(t) = iE{t)e"^^*T - lE* {t)e-"i^'r 



(6) 



The non-equilibrium dynamics is contained in the factor 
E{t)E*{t'), while 



^-W = ([rn<),r(o)]; 



(7) 



becomes invariant under time translation: now the aver- 
age (...) = Trpo(- • ■) is taken with respect to the unper- 
turbed equilibrium density matrix pg cx e^'^^" (all our 
results are valid at arbitrary T = We denote by 

X^(t) = 0{t)X^{t), such that ReX'^iuj) = X_(w). Let 
us first consider the case when V{t) = and \'y{t)\ ~ 1, 
thus E(ui) = 2Tr6{uj). Then, formally, we have: I{qV) — 
—qX_{—qV)/TT. Notice that / is measured as a function 
of V, but we keep q as V appears always multiplied by q, 
when compared to frequency [see Eq. ([8| for example] . It is 
easy to show, by a spectral decomposition, that /(O) = 0, 
which conforms with our assumption that we ignore the 
super-current. However, / is not necessarily an odd func- 
tion of V. We can as well show that the zero-frequency 
noise is given by Eq. ([T]) universally. The finite-frequency 
(FF) noise will be expressed in Ref. flS'. 

Out-of-equilibrium AC conductance: G{qV;uj) 



25 can be expressed either through an exact 
formula [26j 27 adapted for tunneling |13 , or 
from a perturbative computation: 2'KU)G{qV; uj) = 
g2 [X^{^qV + uj)- X^* {-qV - uj) - 2ilm X'^i-qV)] . 
An interesting consequence is that its dissipative part 



ReG{qV;uj) 



_q_ 

2uj 



[I{qV + Uj)-I{qV -u)]. 



(8) 



can be expressed solely in terms of the DC current. This 
is obviously satisfied in the Ohmic regime, but becomes 
less trivial in the nonlinear regime. Similar expression, 
restricted by g = e, has been derived in 21 28 . Note, 



the expressions ([T]) and ([s]) are the sort of the fluctuation- 
dissipation relations, being derived perturbatively with re- 
spect to tunneling. Indeed, Re G yields the asymmetric 
part of the FF noise ^[2j2Q^2!T^ , whose full expression can 
be expressed also in terms of the DC current 13 21 . Re- 
call that Eq.Q is valid at arbitrary T, V, w. We will now 
propose three regimes where it can be used to measure q. 

The form ([s]) offers a first alternative: q can be deter- 
mined through the pref actor and simultaneously by the 
structure of the r.h.s. at w = ±qV . For instance, if / shows 
up a ZBA, Re G{qV] w) gets similar anomalies translated 
a.iV — -izuj/q. This method requires that one of the trans- 
lated voltages y ± oj/g is in the domain of non-linearities 
of /, and to reach frequencies uj ~V. 

Experimentally, low lo are easier to access, which yields 
the second regime, uj <C qV , for which Eq. (Isj) reduces to: 



ReG((7F;a;) - G{qV) 



u? d^GjqV) 



(9) 



where G{qV) = G{qV;uj = 0) denotes the DC differential 
conductance. Such a method offers an interesting and di- 
rect access to q, the r.h.s. can be determined from G{qV). 
Contrary to the previous one, it does not necessitate a 
ZBA. Even more, it is valid in a more accessible regime 
where heating can be reduced: qV, oj <^ T. In this case, 
one needs to access weak non-linearities, the second deriva- 
tive of the DC conductance on the r.h.s. of Eq. ^ being 
now taken at ^ = 0. This is experimentally very promis- 
ing and advantageous compared to the shot noise result, 
Eq. ([!]), which requires at least qV T 14 . 

A third regime corresponds to low voltage with respect 
to frequency, qV ^ w, without specifying T. This is a 
linear regime, where the AC conductance depends only on 
UJ. In particular, if / is odd with respect to voltage, Eq. 
(Isl) reduces simply to 



ReG(0;a;) ~ -I{qV = uj). 



(10) 



In order to get a non-trivial dependence on frequency, uj/q 
has to be inside the domain of nonlinearity of /. Then q 
can be found as a scaling factor by comparing the linear 
AC conductance Re G(0; uj) with the nonlinear DC current 
plotted versus V — uj/q and divided by uj/q. 

Photo-assisted Current (PAC): Let us now consider 
the case of a time-dependent "effective" voltage V{t). It 
may arise either from a time-dependent gate voltage, bar- 
rier height, voltage bias, or from all together. We fo- 
cus, for simplicity, on a periodic variation with a unique 



period 2tt/D,, thus f(t) — J2 fn 



where f{t) stands 



for |7(t)|, V{t) or E{t), and ujn = nVt. Multiple peri- 
ods or more general time-dependence will be treated else- 
where fTsl. Before expressing the PAC, let us specify E^ 
in some typical situations. First, letting |7(t)| — 1 and 
V{t) = 2Vn cos(m), thus K = Vn5n-i for n > 0, 



E„ 



Jr. 



f2qVn 



(11) 
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where J„ are Bessel functions. Second, for \j{t)\ = 1 and 
sharp pulses V{t) = 0o J2n 5[t~2im/VL), we obtain: En — 
sin(#o/2)/(7rn - g0o/2). 

We could also consider a situation where the tunneling 
amplitude is periodic in time while the applied voltage is 
constant. It is possible, as in edge states for instance, to 
modulate the gate voltage: Vg{t) = Vg 
weak Vg, one can expand: 7(t) = 1 + Vg[ze 
where z and z' are complex coefficients, thus 



ac modulation is added, V(t) — woe 



■ v-ne 



iflt 



Then 



Vg COS fit 



For 



En = Sn + Vg{zSn-l + z'Sn+l)- 



(12) 



Let us now compute the time-dependent average cur- 
rent, /(t; qV), whose dependence on all En is not specified 
to avoid confusion. Here we focus on its DC component: 
Io{qV) — (ri/27r) J^^^^ I{qV;t)dt, the easiest to measure 
experimentally. It can be written as: 



(13) 



Thus the PAC is a superposition of the DC currents at 
translated voltages V — nQ./q, due to the absorption of 
n photons, with a probability for each process given by 
|i?„p. This generalizes a main result of the PAT 18 21 
obtained for the special case: for g = e, a tunneling am- 
plitude constant in time, and an AC voltage, so that £"„ 
are given by the Bessel functions, Eq. (111. We propose 



two alternatives to use Eq.( 13 1 for the determination of q 



The first one is suited for the case of a ZBA, taking: 

d'^hiqV) = ^|S„|2ayG(g\/-c.„). (14) 



If dvG is peaked around zero, this expression has peaks 
spaced by = O/g. Here T needs to be small enough 
compared to ri/g, such that it does not smear them. 

The second alternative does not require any ZBA, nei- 
ther low temperatures, thus is less restrictive. It is ob- 
tained by taking the differential with respect to the AC 
components of V, using Eq. (16): 



SVnSV^n 



-2IoiqV) + loiqV + UJn) + h{qV - W„). 

(15) 

Here we have used a general exact property of En (assum- 
ing \"f{t)\ does not depend on V): 



SEn 
SVn,. 



— En 



(16) 



Eq.(15) can be also used to generate higher-order deriva- 



tives. It is powerful as it is valid still at finite and arbitrary 
V. It has not been derived so far within the PAT studies, 
even for noninteracting electrons. We see that once Jq is 
measured as a function of V for a given w„ , comparison of 
both sides allow to determine q. 

Nevertheless this method can be simplified further as 
follows. One starts from a DC voltage V to which a small 



one measures the PAC Io{qV), and takes its second differ- 
ential with respect to vn and w_n- According to Eq. (15) 



in the limit of vanishing V the currents appearing on the 
r.h.s. are nothing but DC currents: 



n2 



S'loiqV) 



SvnSvc 



-2i{qv) + i{qv + + i{qv - n), 



(17) 

A similar formula was derived for q = e in Ref . [2l] . The 
r.h.s. can be obtained from the DC current by shifting the 
argument V by ±f2/q. Thus, for a given frequency the 
quasi-particle charge g is a scaling parameter which makes 
both functions of V to coincide. For this, at least of the 
three voltages, V, V + n/q, and V — il/q, has to be in the 
domain where the DC current is nonlinear. 

Application to the FQHE and quantum wires 
Two edge states in the FQHE are described by a total 
Hamiltonian "Ho in Eq. ([2]), while tunneling between them 
can be induced by pinching off a quantum point contact 
(QPC). We do not require Ho to be decoupled into two 
Hamiltonians, but only to conserve the charges of both 
edges, Qi and Q2- This is a strength of our approach: even 
Coulomb interactions between the edges, usually present 
at least within the QPC, can be incorporated in Hq. An- 
other strength is that the knowledge of Ho is not required 
in order to test the charge in the bulk, or the effective one, 
contrary to methods using the explicit expressions for the 
current or the noise. Moreover, the charge q which appears 
in Eq.([4| is robust against the coupling to an ohmic envi- 
ronment, contrary to the DC conductance without tunnel- 
ing or the exponents controlling the power law behavior 
(of the current and its fiuctuations) [39) . 

We have shown that different regimes can be investi- 
gated, depending on energy scales and technical limita- 
tions. An additional strength, in comparison to the mea- 
surement of the noise in Eq. ([IJ which requires qV 3> T 
14 , is that we have proposed three methods in the regime 
qV T, which reduces heating. Quite often, the nonlin- 
earity domain of I{qV) corresponds rather to qV ^ T. 

Let us focus on these proposals. Among them, two are 
based on the AC conductance; either it is measured at low 
frequencies compared to T, where one needs to access weak 
non-linearities of the DC conductance, see Eq. ([9|, or, ff 
frequencies higher than qV can be reached, one can use di- 
rectly Eq. (10 1. The third one, based on the PAC, can ex- 
ploit Eq. (|17[), which requires here that V +\^\/q:$> T. If 



one needs in addition T 3> qV to reduce heating, one must 
fulfill \n\ ^ T. This is easily accessible as J7 can be varied 
from the RF up to the optical spectrum, thus can be made 
easily larger then many kelvins. Notice that in this regime, 
the r.h.s. of Eq. ([l7| simplifies to -~2qVG + + I{-n). 



Besides being a zero-frequency measurement, this method 
has another advantage compared to the AC conductance: 
one can control the time-dependence of the tunneling am- 
plitude, through the gate voltage, while keeping only a DC 
applied voltage, yielding for instance the amplitudes En 
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in Eq. (12 1. This has less drawbacks compared to an AC 



voltage which needs to propagate along the edges. 

We have assumed that only one process with a given 
charge q dominates at once. It turns out that the value of 
q depends on the regime one considers. For simple frac- 
tions, v = l/(2m + 1), the tunneling charge is q = ve (or 
q — e) at large (respectively, low) energies compared to 
a non-universal energy scale Tb- One can view this as a 
bunching of 2m -I- 1 quasi-particles below Tg . For other se- 
ries of filling factors, the question is still controversial. In 
particular, Heiblum's group claimed (9j[T5] that bunching 
of the fundamental quasi-particles occurs at low energies 
[22] . For instance, at v — 5/2 (or v = 2/3) one measures 
g = e/4 (respectively, q — e/3) at high energy, while one 
observes rather q = e/2 (respectively, q — 2e/3) at low 
energies. Such bunching was though excluded by Radu et 
al [S]. Indeed, Heiblum's group used the expression for 
noise of non-interacting particles with a charge q, which 
has to be taken with caution. Our proposals, which do 
not require any explicit expressions for the current neither 
noise, would be extremely useful to solve these controver- 
sies, and to show with a clear evidence the shift between 
different charges. 

At resonance, i.e., when the gate voltage is adjusted to 
suppress the tunneling amplitude for one quasi-particle, it 
was predicted that for 1/ — l/(2n -I- 1) tunneling of two 
quasi-particles should dominate, thus doubling q ^ 2q 
10 . This was not observed in experiments, q = ve being 
measured 11 , which was explained in Ref. [sT: still one- 



quasi-particle tunneling dominates. 

Another setup concerns junctions between edges at dif- 
ferent filling factors vi and 1^21 for which it was predicted 



32 33 , for simple fractions, that the junction behaves as 
an edge with effective "filling" v — 2viV2/{yi + 1^2)- If 
one couples a Fermi liquid to an edge state [34], one sets 
formally z/i = 1. This has motivated other theoretical 

35 and experimental investigations [36]. The DC con- 



ductance, expected to be given by e^v/h, is nevertheless 
difficult to observe as it is sensitive to the geometry, deple- 
tion, backscattering, etc. One rather expects the effective 
charge ve to be more robust, thus our methods to be of 
great relevance. Yet another setup consists in coupling an 
edge state to a superconductor [37|38] , where the effective 
charge is given hy q — 2ve for simple fractions [33] . 

Our results apply to quantum wires (QWs) or carbon 
nanotubes (CNTs). It was predicted that an injected elec- 
tron splits into two solitonic charges [33j40j4l]. For tun- 
neling junctions as well as weak barriers, the charge q ap- 
pearing in Eq. ([T]) was expected to be respectively given by 
e and Ke [10| , where if is a parameter depending on inter- 
actions, and is roughly the analog of v. However, contrary 
to the FQHE, the coupling to the reservoirs ^(33] [42] 
changes q — Ke to e in Eq. ([T]) [43] , unless high-frequency 



noise is considered 44 . Even though measuring q — e 
looks trivial, it would be an interesting check of the under- 
lying model, adopted so far by the community. More gen- 
erally, provided that Ho is quadratic in the bosonic fields, 



describing arbitrary range, form and inhomogeneities of 
the interactions and parameters, we can show that the 
weak backscattering term in Eq. Q contains q — Go, the 
dimensionless DC conductance of a system described by 
only 45 . This could be relevant for edge states with 
a quadratic action, for instance if the filling factor in the 
QPC is different from that in the outer edges. 

There is a limitation on lo which matters only for the 
methods based on the AC conductance and for 1-D sys- 
tems. The current operator in the electrodes is, strictly 
speaking, different from the tunneling current operator in 
Eq.(|6|. When its average is measured at some distance 
i, the AC average values of both operators coincide only 
at low frequencies, w ^ wl = v/L where v is the plas- 
mon velocity, otherwise the propagation from the tun- 
neling junction to the contacts have to be incorporated 



26 30 



This condition is easily fulfilled for typical QWs 
or CNTs whose length is of the order of 10/im or less, 
thus TTWi, ~ lOOGHz-lTHz depending on v, while uj can- 
not exceed 20GHz. It is as well the case in the FQHE, 
where for L ~ 5/im, and v ~ 5.10^m/s, ujl — lOQGHz. 
One does not expect the limitation on lu to hold in higher 
dimensions, where the plasmon spectrum typically has a 
large gap. Notice also that QWs and CNTs are examples 
where the non-linearity domain of I is not controlled only 
by T. At V <^ Wi/e, even if T ^ eV, the noninteracting 



26 33 



leads impose their "linearity" 

Conclusion: We have considered a non-linear tunnel- 
ing junction between two arbitrary systems in any dimen- 
sion, with any interactions or disorder, excluding only a 
super-current. We have proposed many methods to access 
the charge q of carriers. The main advantages of our pro- 
posals is that the underlying model does not have to be 
specified, that they are based on a current measurement, 
and that many of them are applicable at low frequency 
or/and voltage. They are based on two quantities, each 
of them can be exploited into two different ways: the AC 
conductance 



25 



and the PAC when either the voltage 
or the tunneling barriers, or even both, have a periodic 
time-dependence. 

On one hand, the methods based on the AC conduc- 
tance are proposed in three situations. The first requires 
a ZBA, LU of the order of qV, and at least one ot V ztuj/q 
in the nonlinearity domain of /. Singularities of Eq. ([s]) at 
w = ±qV gives access to the charge. This feature has been 
observed in the FF noise for simple fractions in the FQHE 
16 29 , 30 . Here we have established its universal exis- 



tence in the dissipative AC conductance, which yields the 
asymmetric part of the FF noise. This is related also to 



the universal appearance of a similar feature in Eq. ( 13 1 



The second regime does not require any ZBA, and holds 
even at to, qV ^ T, the easiest regime to attain exper- 
imentally, which reduces heating effects too. The weak 
non-linearities of G{qV) in Eq. ([9| give access to q. 

The third regime is linear with respect to V , as it cor- 
responds to oj qV, where G{qV;uj) ~ G{0;uj). The 
charge q can be deduced easily from Eq. (10). One needs 
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to vary uj/q in the domain where the DC current is non- 
hnear, which depends on the system at hands. In the 
FQHE, it corresponds to w/g » T, accessible experimen- 
tally as well. Typically, a voltage of iO^eV corresponds 
to ~ lOGHz, and a temperature of 50mK to ~ IGHz. 
Technically, it has become possible to reach a; ~ 20Ghz, 
thus greater than both the typical V and T. 

On the other hand, let us recall the methods based 



on the PAC, expressed through the general formula (13). 



Here corresponds to the period 2'K/fl of the applied 
voltage, and can be made much greater than uj above, at 
which the current is measured. The first method requires 



a ZBA and T <C qV such that the spikes in Eq. (14) are 
not smeared out. As noticed before, this is related to the 
similar feature in the AC conductance, Eq. ([s]), and in 
the FF noise [l^. Nevertheless, Iq is much easier to mea- 
sure compared to the FF noise, being in addition at zero 
frequency. The second method based on the PAC, using 



Eq. (17), does not necessitate a ZBA, but an arbitrary 
non-linear behavior. Even more, it does not require low 
temperature compared to both Q and V, but that only 
one of V,V ± fl/q to be in the nonlinearity domain of /. 

Non-linearity is indeed a crucial ingredient in all our 
methods, which are particularly promising for the FQHE. 
They could be useful for super-conducting structures as 
well; the noise in Eq. ([T]) being very weak, more compli- 
cated setups were adopted to measure q = 2e [2]. Alt er- 



natives for linear systems will be proposed elsewhere 13 



We will also express both the current average and its fluc- 
tuations for any time dependence of the voltage and the 
tunneling amplitude in linear or nonlinear systems. 
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